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Abstract 

We study the conversion efficiency of heteronuclear Feshbach molecules in population imbalanced 
atomic gases formed by ramping the magnetic field adiabatically. We extend the recent work [J. 
E. Williams et al., New J. Phys., 8, 150 (2006)] on the theory of Feshbach molecule formations 
to various combinations of quantum statistics of each atomic component. A simple calculation for 
a harmonically trapped ideal gas is in good agreement with the recent experiment [S. B. Papp 
and C. E. Wieman, Phys. Rev. Lett., 97, 180404 (2006)] without any fitting parameters. We 
also give the conversion efficiency as an explicit function of initial peak phase space density of the 
majority species for population imbalanced gases. In the low-density region where Bose-Einstein 
condensation does not appear, the conversion efficiency is a monotonic function of the initial peak 
phase space density, but independent of statistics of a minority component. The quantum statistics 
of majority atoms has a significant effect on the conversion efficiency. In addition, Bose-Einstein 
condensation of an atomic component is the key element determining the maximum conversion 
efficiency. 

PACS numbers: 03.75.Nt 
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I. INTRODUCTION 



There have been a number of efforts to produce molecules in the ultracold atomic gases. 
Making use of the Feshbach resonance, homonuclear molecular conversion efficiencies in 
atomic gases of bosonic *^Rb atoms and fermionic atoms were investigated by Hodby et 
al. They found that the conversion efficiencies do not reach 100% even for an adiabatic 
sweep of the magnetic field, but depend on the peak phase space density of initially prepared 
pure atomic gases. 

The ultracold diatomic heteronudear molecule is now drawing attention. For example, 
there is an expectation of realizing quantum computation by making use of the electric 



dipole moments of heteronudear molecules 



2|]. So far several basic experiments have been 



performed in order to realize this idea. Heteronudear Feshbach resonances have been inves- 
tigated in collisions between bosonic ^^Na atoms and fermionic ^Li atoms j^, bosonic ^''Rb 
atoms and fermionic atoms and bosonic ^''Rb atoms and bosonic ^^Rb atoms ^. 
The adiabatic conversion efficiency of heteronudear Feshbach molecules composed of ^''Rb 
and ®^Rb in population imbalanced gas have been investigated in detail j^. Fermionic he ;- 
eronuclear Feshbach molecules via rf association are investigated in a 3D optical lattice joj. 

In order to analyze the experimental data for the molecular conversion efficiency, Hodby 
et al. [V\ devised a stochastic phase space sampling (SPSS) model based on the assumption 
that only a pair of atoms sufficiently close in phase space can form a molecule through 
the adiabatic process. With a single fitting parameter 7 (a cutoff-radius in phase space), 
which determines how close two atoms should be in the phase space to form a molecule, 
the phenomenological SPSS model predicted that the molecular conversion efficiency is a 
universal function of the initial peak phase space density. The prediction from the SPSS 
model agreed notably well with the experimental data [ll]. The SPSS model was also used to 



analyze the conversion efficiencies ^''Rb-^^Rb Feshbach molecules [5[ , which was in agreement 
with the experimental data. 

On the other hand, Williams et al. [t^ provided the theory of Feshbach molecule forma- 
tion based on the coupled atom-molecule Boltzmann equation. They derived from the first 
principles the universal functional form for the conversion efficiency without any fitting pa- 
rameters. This approach is appropriate for the regime tramp ^ Tcoh where multiple collisions 
occur. Here, tramp is an upper bound on the ramp time and Tcoi is the average time between 



2 



collisions. 

The coupled atom-molecule Boltzmann approach provided a result that the exchange of 
atoms and molecules ceases when the conservation laws of momentum and energy via the 
atom- molecule conversion cannot be satisfied {3, Q]- When we neglect the energy shift by the 
real part of the self-energy, these conservation laws are not satisfied for negative detunings. 
The detuning is the difference between an energy of the atomic dissociation continuum and 
the lowest single particle energy of the molecule, which is controlled by the magnetic field 
due to the Zeeman shift. 

n 

According to Williams et al. [7], the conversion mechanism can be described as follows. 
When the starting point is the positive detuning far from the resonance with pure atomic 
gas, Feshbach molecules are formed and decay as the detuning is ramped in the region 6 > 0. 
When the detuning crosses zero and becomes negative, the molecule formation effectively 
halts, which gives rise to the saturation of molecule production because the momentum and 
energy conservation laws are not satisfied. If the density is high and the 3-body recombi- 
nation to form the molecule is effective, the saturation of molecule population would not 
appear, and the molecule formation would take place even in negative detunings. However, 



the saturation of the mo 
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e population against the magnetic field is observed in many 



12, 



13| . Therefore, one is allowed to assume that the sweep 



speed of the magnetic field is fast with respect to the time scale for 3-body collisions and 
only the atom-molecule resonant collisions play a role in producing the molecule in these 
experiments. 

On the other hand, the saturation against the inverse sweep rate is reported in experi- 
ments , [l^ . The saturation of molecular population with slow sweep means that the 
detailed balance is realized. Using the kinetic theory, Williams et al. j3] signified that the 
atom-molecule resonant collisions alone can give rise to thermal equilibrium as well as chem- 

n 

ical equilibrium [8| . Considering the assumption that a sweep is slow with respect to 2-body 
process but it is fast with respect to the time scale for 3-body recombinations, they concluded 
that the detailed balance is only realized for positive detunings only through atom-molecule 
resonant collisions. We note that the production rate of entropy vanishes as long as the slow 
sweep of the magnetic field keeps the system in local equilibrium. Therefore, when the sweep 
is so slow that the molecular conversion rate gets saturated, the process may be regarded 
as adiabatic for positive detunings. For negative detunings, the chemical equilibrium is no 
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longer maintained. 

As a result, the conversion efficiency with adiabatic sweep is determined by the molec- 
ular fraction at zero detuning, which is a function of the physical quantity of initial state 
connected through the equal entropy. This mechanism does not require any fitting param- 
eters. According to this mechanism, the conversion efficiencies plotted against the initial 
peak phase space density [3] agreed well qualitatively with the experiment even though 
using an ideal gas mixture model. Although these saturation mechanisms are derived above 
the superfluid transition temperature, these mechanisms are based on the momentum and 
energy conversion laws so that they are assumed to be quite general. The principle given in 
Ref. 31 extended to the superfiuid phase of a Fermi gas below the transition tempera- 
ture pj;]. Conversion efficiency against the initial temperature [l^ including the resonant 
interaction agreed well qualitatively with the experiment 



Within the classical gas approximation, the explicit formula 



was given as a function of the initial peak phase space density [t], which shows a qualitative 
agreement with the experiment consistent with the SPSS model. 

In this paper, we study the conversion efficiency of the heteronuclear Feshbach molecule 
in population imbalanced gases, extending the theoretical model provided by Williams et 



"or the conversion efficiency 



J. 



al. fTll. We will show that our calculation agrees with the experimental result by Papp and 
Wieman [sl without any fitting parameters. An exquisite point of this theoretical model is 
that one does not need any fitting parameters so that one can make quantitative prediction 
for the conversion efficiencies which have not been measured yet. Thus, we investigate 
various combinations of quantum statistics of atoms. All combinations studied in this paper 
are collected on the table below. Symbols "B" and "F" represent the quantum statistics of 
boson and fermion, and symbols ">", "<" and "m" represent the majority atomic component, 
the minority atomic component, and the molecular component, 

'(1) B> + B<^B^, 
(2) B> + F<^F, 



(3) F> + B<^F, 



^ (4) F> + F< ^ B^. 

We also give the conversion efficiency as an explicit function of the initial peak phase space 
density of the majority atomic component. We will show that in the low-density region 
where Bose-Einstein condensation does not appear the conversion efficiency is a monotonic 



function of the initial peak phase space density, but independent of statistics of the minority 
component. 

For homonuclear Feshbach molecules, the trend of conversion efficiencies calculated by 



Williams et al. 
in Refs. jl] and 



was qualitatively consistent with the SPSS model within the rarige plotted 
Tj. However, we will show that the model by Williams et al. ^ leads to 
completely different results from the phenomenological SPSS model at low temperatures, 
when applied to heteronuclear Feshbach molecules. In the case of bosonic homonuclear 
Feshbach molecules composed of bosonic atoms, the trend of conversion efficiency is also 
expected to have a completely different result from the SPSS model in the presence of Bose- 
Einstein condensation. Although the SPSS model has shown agreement with the experiments 
so far, the problem with this model with regard to heteronuclear Feshbach molecules is that 
the atoms of different species are disposed in the same phase space, although they should 
be essentially in own phase spaces. 

II. EQUILIBRIUM THEORY 

In this Section, we present the equilibrium theory for the ideal gas mixture composed of 
two species of atoms and heteronuclear Feshbach molecules, assuming that atomic popula- 
tions are imbalanced. 

We denote the numbers of majority atoms, minority atoms and molecules as A^>, A^< and 
A^m, respectively. The following two constraints should be satisfied in these systems. The 
first constraint is the conservation of the total number of particles: 

iVtot = Ar> + Ar< + 2Ar^. (2) 

The second constraint is the conservation of the ratio: 

where a is determined by initial ratio of A^<,ini and A^>,ini (where there is no molecular 
component), i.e., iV< ini/A^>,ini = By definition of the majority and minority components, 
one has a < 1. The second constraint can also be written as 

A^< - aN> + (1 - a)Nra = 0. (4) 
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In order to impose two constraints in Eqs. ([2]) and (jl]) described above, we deal with 
the grand canonical Hamiltonian with introducing two Lagrange multipliers fii and fi2- The 
partition function S for this grand canonical Hamiltonian is defined by 



Tr 



exp {^0 - /ii (iV> + iV< + 2iV^) - ^2 [iV< - aiV> + (1 - a)iV^] })] , (5) 



where = l/k^T with /cb and T being the Boltzmann's constant and temperature. The 
Hamiltonian Hq for the ideal gas mixtures is given by 

i j k 

where ef, ef and are the single-particle energies, and hf , hj, and are the number 
operators. The lowest energy of the molecule is displaced relative to that of the atomic dis- 
sociation continuum by the detuning 6, which is controlled by the magnetic field. Operators 
of atomic and molecular populations are given by 

i j k 

The thermodynamic potential VL is defined hy Vt = — fceTlnS. We derive two equations 
representing two constraints from the thermodynamic potential Vt. The total number of 
particles is given by A^tot = —dQ/dfii, which leads to 



The equation representing the second constraint is given by = —dQ/dfi2, which leads to 
= V i "V — r-h + (1 - «) V — r^rs • (9) 

In Eqs. dH]) and iQ, the negative sign is for bosons and the positive sign is for fermions. The 
fugacities of the majority atoms, the minority atoms and the heteronuclear molecules are 
defined by 

These fugacities are related to each other through 

z^z^ = z^e'l^^^. (11) 



15|. 



The relation at zero detuning S = given by -2>-2< = will play an important role in 
considering the molecular conversion efficiency. 

We consider a gas confined in an anisotropic harmonic trap with frequencies {cj* , Uy, ul}, 
where i represents the component indexes {>, <, m}. For later use, we define the ratio of the 
trap frequencies 7< = COy/Cd^ and 7m = uj^/uj^, where cUi = (cu^cu^cj*)-^/^. In optical dipo 
traps, the trap frequencies are determined by the polarizability and the particle masses 

As in Ref. [l^ , we assume ksT ^ huji and thus use the semi-classical formulation, which 
replaces the sum over discrete states in Eqs. ([8]) and Q) by an integral. The density of states 
Pi{e) in a harmonic trap is given by Pi{e) = e"^ /[2{}TWi)^]. The noncondensed population Nb 
for bosons and the population N-p for fermions are represented by the Bose integral Qn{z) 
and the Fermi integral J-'n{z) respectively, which are given by 



M^^^^jg,^.^, (12) 

where z is a fugacity, and Co is the geometric average of the trap frequency u = {u^UyUzy^^. 
The Bose and Fermi integrals are defined by 

1 poo n-lJ 1 poo n-lJ 

^n{z) ^ — - / -T— -, (13) 



' T{n) Jo z~^e^-V ' ^{n) Jo z-^e"^ + 1 

where T{n) is the gamma function. The Bose function is equal to the Riemann-zeta function 
({n) when z = 1, i.e., Qn{z = 1) = (i^^)- 

The entropy of noncondensed bosons and the entropy of fermions Sp are given by [l^ 



Sb(z)= fceiVB 



4^-lnz 



Sv(z)= k^Np 



4 ) [ -\nz 
J's z) 



(14) 



In the low temperature limit, they are approximated by 



,„.«e(|^)'c(4,, S...v(f)"|^. (15) 

We assume that the initially prepared state has a large and positive detuning so that no 
molecules exist in the initial state. The initial populations of atoms satisfying two constraints 
in Eqs. ([2]) and (jlj) are given by 

1 a 
N>,ini = -—Ntot, Ar<,im = ^— A^tot. (16) 
1 + a 1 + a 

According to the model for the molecular conversion efficiency proposed in Ref. [?], we 
define the molecular conversion efficiency Xo as a fraction of a molecular population at 



zero detuning 5 = 0, given by 



Xo 



N^{6 = 0) iV^(5 = 0) 1 + a 



N. 



<,mi 



tot 



a 



(17) 



III. BOSONIC HETERONUCLEAR MOLECULES COMPOSED OF BOSONIC 
ATOMS 



In this Section, we consider the first case: {B> + B< ^ Bm}. To make contact with 



the experiment of Ref. [5j, we relabel two atomic components ^"^Rb and ^^Rb as 87 and 85 
respectively, where ^^Rb is a majority component and ^^Rb is a minority component. The 
value of a is chosen from experimental data; we^ adopt a = 2/15, assuming A^85,ini = 40, 000 
and NsT^ini = 300, 000 reported by Papp et al. [5] to be typically produced. 
We note that the two constraints in Eqs. ([2]) and (jlj) are given as 



tot 



knT 



•87 



^3(^87) 



+ 



iVf + 



+2 



knT 



fkU; 



•85 



^3(^85) 



^3(2m) 



iVf + 



'85 



^3(^85) 



a 



iVf + 



kuT 



87 



^3(^87) 



+ (1 



a 



, kjiT 

fVjJm 



(18) 



where A*"^^, A*"^^ and N™ are condensed populations of the respective components. 

According to the experiment Q|, ^^Rb is a majority component, compared with ^^Rb and 
the heteronuclear molecule. This makes the Bose-Einstein condensation (BEC) transition 
temperature of ^''Rb atoms, which we define as T^'^ , higher than that of ®^Rb atoms and 
heteronuclear molecules. The condition for BEC of ^''Rb component is given by 



(19) 



We note that the fugacities ^§5 and given in Eq. (ITUl) are always less than 1, because the 
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distribution functions must be positive. Below T^'^ , these conditions reduce to 

( {l + a)fi2 < 0, 
1 (l + a)/i2-5 < 0, 



(20) 



where we have used the condition /ii = a/i2. 

We shall determine the phase diagram in 6-T plane. Below T^f^, there are some possible 
cases where ^^Rb atoms or heteronuclear molecules, or both could be condensed. We consider 
three separate cases; 6 > 0, 6 < and 6 = 0, because populations of bosonic atoms ^''Rb 
and ^^Rb and bosonic heteronuclear molecules are sensitive to the sign of the detuning 6. 

We first consider the case S > 0. According to the inequalities in Eq. (!20|) and 5 > 0, we 
find that 6 — {1 + a)fi2 cannot reach 0. In other words, heteronuclear molecule can never 
be condensed {N^ = 0). On the other hand, BEG of ^^Rb atoms takes place when fi2 = 0. 
This is reasonable because the lowest energy of heteronuclear molecules is higher than that 
of *^Rb atoms. We determine below which BEG of *^Rb atoms takes place from the 
following equation. 



with Zm = exp[— 5/(A;bT+)]. Below T+, one has the following exact forms for two condensate 
populations. 



where 2^ = exp [— 5/(A;bT)]. We note that in the limit T 0, the above condensate 
populations correspond to the respective initial populations so that one has no heteronuclear 
molecules at 5 > in the limit T ^ 0. This is specific to the case of a Bose-Bose ideal gas 
mixture, where the ground state of the system only depends on the sign of the detuning. In 
contrast, the ground state of a gas including a fermionic component, which we will discuss in 
the following sections, is dominated by Fermi energy. One can have heteronuclear molecules 
even in the limit T ^ at positive detunings in such systems. Strictly speaking, these 
results in the limit T ^ are due to our treatment assuming k^T ^ hwi. In the limit 
T — > satisfying k^T < hui, the quantized discrete energy spectrum of each component 
cannot be ignored. 




(21) 




(22) 
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Secondly, we consider the case S < 0. According to the inequahties in Eq. (120|) and S < 0, 
we find that /i2 cannot reach 0. In other words, ^^Rb atoms cannot be condensed (A^^^ = 0), 
while BEC of heteronuclear molecules can take place under the condition 6 = {1 + a)fi2- 
This is reasonable because the lowest energy of heteronuclear molecules is lower than that 
of ^^Rb atoms. One can determine T~ below which BEC of heteronuclear molecules takes 
place by solving the following equation, 

j-^j [7mC(3) + 7156^3(^85)] -aArtot = 0, (23) 
where 2^5 = exp [6/kBT~]. Below , we obtain the exact forms of two condensate popula- 



tions. 



with = exp [5/{kBT)]. 

Finally, we consider the case 5 = 0. Although BEC of ^^Rb atoms and heteronuclear 
molecules can take place below T^''' when /ii = /i2 = 5 = is satisfied, this case has a problem 
in that condensate fractions cannot be uniquely determined because of the degeneracy of 
^^Rb atoms, ^^Rb atoms and heteronuclear molecules. The transition temperature T° below 
which BEC of ^^Rb atoms and heteronuclear molecules appears, with ^''Rb atoms being 
condensed, is given by 

The condensate population N^"^ at T = T° is given by 

Nr = iVtot - (^) ' C(3) (1 + 78^ + 27^!.) . (26) 

So far we have not specified trap frequency for each component. A trap frequency depends 
on a mass and a polarizability of a particle. The polarizability of a molecule is expressed by 
the sum of the atomic polarizabilities, when the molecular internuclear separation is large 



(25) 



compared to the mean atomic radii ISj. Porlarizabilities of ^^Rb atoms and ^^Rb atoms 
are equal, because atoms are isotopic. In this isotopic case, ratios of trap frequencies 735 
and 7m are expressed by masses of two components according to simple formulae given in 

10 
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FIG. 1: The phase diagram of Rb atoms, Rb atoms and heteronuclear molecules. We assume 
the ratio a = 2/15, and equal trap frequencies (Ds? = wss = '^m- 

Ref. [13] . Two atomic masses are different, but the difference of masses are so small that 
the ratios of trap frequencies are almost unity; 735 ~ 0.99 and 7m ~ 0.99. As a result, one 
is allowed to assume that all trap frequencies are equal cjgy = cjgs = uj^. We will use equal 
trap frequencies henceforth. 

We plot the phase diagram for an ideal gas mixture of *^Rb atoms, '^^Rb atoms and 
heteronuclear molecules in Fig. [1] against the detuning 5 and the temperature T. Fig. [2] 
shows population of each component. Figs. [2] (A), [2] (B) and[2](C) show condensate fractions 
of ^^Rb atoms, ^^Rb atoms and heteronuclear Feshbach molecules, defined by N^'^ /Ntot, 
N^^Ntot and N^/Ntot- Figs. [2] (D), [2] (E) and [2] (F) show fractions of «^Rb atoms, ^^Rb 
atoms and heteronuclear Feshbach molecules, defined by Nsr/Ntot, Ags/A^tot and N^/Ntot- 
As discussed above, populations of ^''Rb atoms, ®^Rb atoms and heteronuclear molecules 
change dramatically and discontinuously on both sides of the 5 = at low temperatures. 

In order to form Feshbach molecules, an adiabatic ramp of the magnetic field is typically 
used. Throughout this process, the entropy should be conserved. The total entropy of this 
system is given by 

S = Sb{zs7) + Sb{zs5) + Ssiz^), (27) 
where the exphcit formula of Sb is given in Eq. (fT^ . In Fig. [3l we plot the contours of 
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FIG. 2: (Color online) (A) The condensate fraction of ^'^Rb atoms: N^"^ /Ntot. (B) The condensate 
fraction of ®^Rb atoms: N^^/N^of (C) The condensate fraction of molecules: N^/Ntot- (D) The 
fraction of ^'^Rb atoms: Nsj/Ntot- (E) The fraction of *^Rb atoms: Ns^/Ntot- (F) The fraction of 
molecules: A^m/-^tot- We assume the ratio a = 2/15, and equal trap frequencies wsy = (Dss = Wm- 
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FIG. 3: (Color online) The contours of constant entropy. We assume the ratio a = 2/15, and 
equal trap frequencies 0)37 = Qs5 = <^m- 

constant entropy assuming the system traverses through the adiabatic ramp. 

The adiabatic process connects an initial state with the final state by the equal entropy. 
Initial numbers of ^^Rb atoms and ^^Rb atoms are given by 

1 



87,ini 



tot ) 



a 



85,ini 



tot- 



We define the initial transition temperatures for ®^Rb atoms and ®^Rb atoms by 



(28) 



'^B-tc.ini 



87 



87,im 



nl/3 



'^B-ic.in 



85 



-^85,im 



1/3 



(29) 



The transition temperatures satisfy T^J^^ > T^f^^ since ^^Rb is a majority component. For 
T'c fni ^ ^ini ^ Li' ohIj BEC of ^''Rb atoms occurs under the condition given by /Ji = a^2- 
On the other hand, both ^^Rb atoms and ^^Rb atoms arc condensed for < Tini < T^\-ai 
under the condition /ii — 112 — 0. In this temperature region, the condensate populations 
are given by 

3 



l + a 
a 



1 + a 



Ntot- 
Ntot- 



^B^ini 
JiLJ. 



'87 



85 



C(3), 

as). 



(30) 
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We now determine the conversion efficiency from the heteronuclear molecule fraction at 
zero detuning. As noted above, in the case of Bose-Bose ideal gas mixtures, one cannot 
uniquely determine the heteronuclear molecule fraction at 5 = below Tc. We thus adopt 
the molecular conversion efficiency as the heteronuclear molecule fraction just above zero 
detuning 5^0+; 

^ N^{5 ^ 0+) ^ N^jS ^ 0+) 1 + g 

iV85,i,i iVtot « ■ ^ ^ 

This treatment is appropriate since the adiabatic sweep starts at 5 > side and a formation 

Q 

rate and a decay rate vanish at 5 = |7|. 

By equating a final entropy S'f (Tf, 5^0+) just above the zero detuning at 5 — 0+ with 
an initial entropy S'ini(Tini) = 5'B(-287,ini) + 'S'B(2;85,ini), we obtain the relation between the final 
temperature Tf and the initial temperature Tin;. In the low temperature limit, we obtain the 
explicit formula, 

1,3 \ 1/3 



Tf = ^ ^ 3 3 Ti,i. (32) 

.1 +785 +7^' 



Fig. H] shows the relation between the initial temperature and the final temperature de- 
termined from Sf(Tf,6 0+) = 5'ini(Tini). The solid line is the numerical result and the 
dotted line is the low temperature limit given in Eq. (1321) . In the high temperature limit, the 
heteronuclear molecular population is so small that the contribution from the heteronuclear 
molecule entropy to the total entropy is small. As a result, one has Tf ^ Tini in the high 
temperature region. 

Fig. [5] shows the molecular conversion efficiency as a function of the initial temperature. 
The dots are the experimental data 5l . The dot-dashed line is the result of our calculation. 

n 

The solid line shows the result of the SPSS model and the dashed lines are its uncertainty [5| , 
where we note that a Maxwell-Boltzmann distribution is always used for ^^Rb atoms in their 
simulation j^. Our calculation agrees quite well with the behavior of the experiment {sl 
without any ffiting parameters. 

In the experiment 5j, the ratio a is changed in order to vary T/T^J^^. However, the 
qualitative behavior of our calculation is unchanged even if we change the value a. Especially, 
the conversion efficiency has a plateau at Xo = '^85/('^85 + '^m); where ^'^Rb atoms are 
condensed. This value does not depend on a. The plateau is also seen in the experimental 
data. 
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FIG. 4: (Color online) The relation between the initial temperature {5 — oo) and the final 
temperature {6 0"*"). The solid line is the numerical result, while the dotted line is the result 
of the low temperature limit given in Eq. (|32p . We assume the ratio a = 2/15 and equal trap 
frequencies (Dgy = Os5 = ^m- 

The origin of the deviation of our result from the experimental result may be due to our 
treatment using the ideal gas mixture. Including the mean-field interaction for the majority 
component ^^Rb atoms in our model will shift the results. 

Below the plateau region, the conversion efficiency starts to decrease when both *^Rb and 
^^Rb atoms become Bose-condensed. Comparing our result with the experimental data, the 
data point at the lowest temperature (xo — 40%) might be considered as the signal of ^^Rb 
BEC. We write down the conversion efficiency xo fhe region where both ^''Rb and ^^Rb 
components are condensed, as a function of the initial temperature Tini, as follows 

= fMiV ^(3)^1±^ = (1 + 7|5)7^ 
° V / iVtot a 1 + 7I5 + 7^ 

For Tini = 0, one cannot convert any atoms into heteronuclear Feshbach molecules. 

Although both conversion efficiencies in our model (dot-dashed line) and in the SPSS 

n 

model [5j (solid line) drop at low temperatures as shown in Fig. [5l we note that mechanisms 
of the decrease in the efficiency are fundamentally different. In our model, molecules are 
not condensed in the region where the efficiency decreases, while ^^Rb atoms are condensed 



A^B^ini 



87 



C(3) 



1 1 + a 



a 



(33) 
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since S ^ 0~^. The total energy becomes lower when heteronuclear molecules decay into 
atoms and become Bose-condensed, because the lowest energy of ^^Rb atoms is lower than 
that of heteronuclear molecules for 6-^0^. This leads to the decrease of the conversion 
efficiency in our model. On the other hand, the decrease of the efficiency in the SPSS model 
is entirely due to the use of the Maxwell-Boltzmann distribution for a minority component 
of *^Rb atoms at all temperatures. In the region where BEG of ^^Rb atoms takes place, 
^^Rb atoms are localized in the phase space, while ^^Rb atoms spread in the phase space 
due to the classical gas distribution. Thus, in the SPSS model, it becomes difficult for ^^Rb 
atoms to find partners with ^''Rb atoms as temperature decreases, because of the principle 
of a monogamy clause for a proximity atomic pair in phase space. This leads to the decrease 
of the conversion efficiency in the SPSS model. 

On the other hand, when ^^Rb gas is Bose condensed, the SPSS model will predict that 
the molecular conversion efficiency quickly approaches 100% as the gases go to the ground 
state [51]. This is because every atom occupies the same state in phase space if the system 
is in the ground state, and ^^Rb atoms can easily find partners with ^^Rb atoms close to 
them in phase space due to the principle of the SPSS model. We note that our theory would 
have also predicted that the conversion efficiency reaches 100% if we defined the conversion 
efficiency as xo = ^m(5 — ^ 0^)/Ns5,im just below zero detuning. However, this treatment is 
inappropriate to the argument based on the energy and momentum conservation proposed 
by Williams et al. [3]. 

The conversion efficiency xo in the plateau region, where only ^^Rb atoms are condensed, 
can be determined as follows. We recall from Eq. (ITTl) that the fugacities satisfy ^87 ^85 = z-^i 
at zero detuning. When ®^Rb atoms are condensed, i.e., z^'j = 1, one has the relation 
2^85 = -2m, which makes the population ratio of ^^Rb atoms and heteronuclear molecules 
independent of the temperature. As a result, the conversion efficiency Xo keeps the constant 
value, exhibiting the plateau. This maximum conversion does not depend on the atomic 
population ratio a, but does depend on the trap frequencies for minority component and 
molecule; 



XO,max = _3 ^^3 - (34) 

^85 + 

Even if we regard ®^Rb atoms and heteronuclear molecules as classical gases, this mechanism 



does not change and the conversion efficiency has a plateau when ^^Rb atoms are condensed. 
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FIG. 5: (Color online) The conversion efficiency as a function of the initial temperature. We 
assume the ratio a = 2/15, and equal trap frequencies (Dg? = cDss = ^m- 

i.e., = 1. This result is also different from the result of the SPSS model, which shows 
that the efficiency decreases in the low temperature region when ®^Rb atoms are regarded 
as the classical gas. From Eq. (IMIl . one sees that the maximum conversion efficiency in our 
model is given by Xo,max = 50%, for equal trap frequencies ljs7 = (^85 = ^m- 

Before closing this section, we comment on the conversion efficiencies of bosonic homonu- 
clear Feshbach molecules composed of bosonic atoms. In this case, the theoretical result 



using the classical gas approximation based on the theoretical model by Williams et al. [?] 
is qualitatively consistent with the result of experiment and the SPSS model [1]. However, 
when we extend Ref. Q to the low temperature regime where Bose-Einsten condensation 
appears, the result is expected to be also completely different from the result by the SPSS 
model. In the SPSS model, one expects a complete conversion as long as Bose-Einstein 
condensation of atomic component appears. On the other hand, applying the theoretical 
model by Williams et al. [7] , one finds that the conversion efficiency does not have a plateau 
but starts to decrease when BEC of atoms appears. As a result, one cannot convert any 
atoms into molecules in the limit Tjni 0. This result can be obtained from the fugacity 
relation at zero detuning = z^, where Za is the fugacity of atomic component, and the 
definition of the molecular conversion efficiency using the molecular population just above 
zero detuning N^{6 — ^ 0+). 
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IV. GAS OF MAJORITY BOSONIC ATOM, MINORITY FERMIONIC ATOM 
AND ITS HETERONUCLEAR MOLECULE 



In this section, we discuss the conversion efficiency of heteronuclear Feshbach molecules 
composed of bosonic and fermionic atoms. We ffist note that the phase diagram of an ideal 
Bose- Fermi mixture gas with a majority bosonic component is not smoothly connected to 
the case of a majority fermionic component. When a majority atomic component is bosonic, 
not all bosonic atoms are converted into heteronuclear molecules and thus the remaining 
atoms can undergo Bose-Einstein condensation at any detuning. In contrast, in the case of 
a majority fermionic component, the transition temperature of Bose-Einstein condensation 
drops to zero at a certain detuning where all bosonic atoms are converted into heteronuclear 
molecules. The phase diagram of a population balanced gas of Bose- Fermi mixture belongs to 
the latter case. Morales et al. studied populations and adiabatic trajectories in 5-T plane for 
an ideal trapped gas of bosonic and fermionic atoms and the heteronuclear molecules formed 
from these atomic components 15^. Their numerical results assumed equal populations in 
the two atomic components and different trap frequencies. 

In this section, we discuss the case with a majority bosonic atom: {B> -|- F< Fm}, 
following the same procedure as in the previous section. As noted above, the calculations 
given in this section is only valid for a < 1. On the other hand, the population balanced 
gas of Bose-Fermi mixture characterized by a = 1 is smoothly connected to the case with a 
majority fermionic component: {F> + B< ^ Fm}, which we will discuss in the next section. 

Under the two constraints in Eqs. ((21) and (HI), the populations of three components are 
given by 

where is the number of condensate atoms in a bosonic majority component. The condi- 
tion for a majority component BEG is yUi = a^2- Fig. [6] shows the phase diagram showing 
the boundary of condensate phase of a bosonic majority component. Fig. [7] shows fraction 
of each component, where we assume a = 2/15. We also assume equal trap frequencies 
ti}> = c<}< = cUjn for simplicity. Figs. [7] (A), [7] (B) and [7] (C) show fractions of the ma- 
jority atomic, the minority atomic and the heteronuclear molecular component, defined by 
A^>/A^tot, N^/Ntot and N^/Ntot- Fig. [7| (D) shows the condensed fraction of the majority 
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FIG. 6: The phase diagram of the transition temperature of a bosonic majority component. We 
assume the ratio a = 2/15, and equal trap frequencies (D> = (D< = (Dm- 



atomic component defined by /Ntot- Fig. [8] shows the contours of constant entropy as- 
suming the system traverses through the adiabatic ramp, where the total entropy is given 
by S* = 5'b(2;>) + Sf{z<) + Sf{z^). 

Under the condition /xi = afj,2, the BEC transition temperature T^. at zero detuning 5 = 
is given by 



7^ + (1 - a)'y 



1 



tot 



1/3 



(36) 



7| + 7;L (1 + «)C(3)- 
Right at zero detuning 5 = and below T^, the populations of the bosonic majority com- 
ponent BEC A^^-^, the fermionic heteronuclear molecule N^, and the fermionic minority 
component A^< are given by 

1 73+ (1-0)73 



N>{6 = 0) 
N^{6 = 0) 



1 + a 
a 



y3 
/ m 



1 -h a 7^ 7, 



7< + 7n 
^3 
!m AT 

^Vtot, 



"^Vtot 



<,ini 



- iVm(5 = 0) 



a 



kBT 



C(3), 



(37) 



tot- 



1 a 73, 73 

In the low temperature limit, one can obtain an analytical expression for the final entropy 
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FIG. 7: (Color online) (A) The fraction of a bosonic majority component: N^/Ntot. (B) The 
fraction of a fermionic minority component: N<^/Ntot- (C) The fraction of a fermionic heteronuclear 
molecular component: N^/Ntof (D) The condensate fraction of a bosonic majority component: 
/Ntot- We assume the ratio a = 2/15, and equal trap frequencies ci)> = (D< 



UJ„ 



at zero detuning 5 = as 



/CB7r'(7< + 7n 



C(4) + ^Bvr^ 



(iV< 



6 
1 



1/3 



6 V 1 + a 7| + 7m 



tot 



1/3 



1/3 



(38) 



where the approximation leading to the last line of Eq. (l38l) is valid when the final temper- 
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FIG. 8: (Color online) The contours of constant entropy. We assume the ratio a = 2/15, and 
equal trap frequencies a)> = a;< = Ujn- 



ature satisfies the condition 



(39) 



In a pure atomic gas before a sweep of the magnetic field, initial populations in normal 
phase are given by 



N- 



1 



> ,ini ^ , ' tot 

1 + a 



^3(2;>,ini), 



^ 3(^<,im)- 



(40) 



1 + q; y tiw<: 

The transition temperature Tc^mi and a condensate population N^^^^^ below Tc^mi are given 
by 

VC(3)j ' 



^B^Cini — ^> I 



N> ■ 

cam 



1 



1 + a 



-7V< 



tot 



^B^ini 



(41) 



C(3). 
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In the low temperature limit, the initial entropy is given by 



5'ini — >S'B(2;>,ini) + 5'f(^<, 



Ah 



If a 



C(4) + k^T:' 

?1 1/3 



6 V 1 + a 



tot 



(A^<,ini)^ 

6 

^B^ini 



1/3 



(42) 



where the approximation leading to the last line of Eq. fj42|) is valid when the initial tem- 
perature satisfies the condition 




1/2 fc- 



(43) 



Connecting the final entropy at zero detuning 5 = with the initial entropy, we obtain 
the relation between an initial temperature Tin; and a final temperature Tf. Especially in 
the low temperature limit, using Eqs. (|38l) and (H2l) . we obtain the explicit relation 



7< 



(^3+^3^)1/3 



T- 



(44) 



In the high temperature limit, the heteronuclear molecular population is so small that the 
contribution of the heteronuclear molecule entropy to the total entropy is very small. As a 
result, one has Tf ^ Tin; in the high temperature region. Fig. [9] shows the numerical result 
for the final temperature as a function of the initial temperature. For comparison, we also 
plot the analytical result in the low temperature limit given in Eq. (jSj). 

Fig. [TOl shows the conversion efficiency as a function of the initial temperature, calculated 
from the molecular population at zero detuning 5 = 0. The molecular conversion efficiency 
Xo ill a plateau region, where the bosonic majority component is Bose condensed, is given 
by 

N^{5 = 0) 1 + a 



r)3 



(45) 



iVtot a 

As shown in Fig. [TUl this value is the maximum conversion efficiency, which does not depend 
on the number ratio a. From the relation between fugacities z^z^ = z^ at zero detuning 
5 = 0, we find z^ = z^ when the majority component is Bose condensed, i.e., z> = 1. This 
makes number ratio of the minority component and the heteronuclear molecule constant, 
and thus the molecular conversion efficiency exhibits a plateau. 
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FIG. 9: (Color online) The relation between the initial temperature and the final temperature at 
6 = 0. The solid line is the numerical result and the dotted line is the analytical result in the low 
temperature limit given in Eq. (j44p . We assume the ratio a = 2/15, and equal trap frequencies 

We note that this result is not expected from the analysis using the SPSS model. As 
the system approaches the ground state, bosonic atoms occupy the same state in the phase 
space. In contrast, fermionic atoms spread forming the Fermi surface due to Pauli exclusion 
principle. At low temperatures where the quantum degeneracy appears, the SPSS model 
predicts that the conversion efficiency decreases as the temperature goes to zero, because the 
number of atoms that can find a pair decreases due to a monogamy clause for a proximity 
atomic pair in phase space. The prescription disposing the atoms of different species in the 
same phase space leads to this result. 



V. GAS OF MAJORITY FERMIONIC ATOM, MINORITY BOSONIC ATOM 
AND ITS HETERONUCLEAR MOLECULE 

In this section, we consider the third case, with a majority component of fermionic atoms: 
{F> + B< Fill}. Under the two constraints in Eqs. ([2]) and (jl]), the populations of three 
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FIG. 10: The molecular conversion efficiency as a function of the initial temperature. We assume 
the ratio a = 2/15, and equal trap frequencies (D> = w< = (Dm- 

components are given by 

N> = (^)'-^3(^>), iV< = N< + (^^y^3(^<), = (^^y ^3(^^), (46) 

where A*"^*^ is the condensate population of minority atoms. The condition for a minority 
component BEC is given by /ii = —fi2- Fig. [TT] shows the phase diagram showing the phase 
boundary for BEC in a bosonic minority component, where we assume a = 3/4 and equal 
trap frequencies = = for simplicity. Fig. [12] shows fraction of each component. 
Figs. [12] (A), [12] (B) and [12] (C) show fractions of the majority atomic, the minority atomic 
and the heteronuclear molecular components defined by Ny/Ntot, N^/Ntot and N^/Ntot- 
Fig. [12] (D) shows the condensed fraction of the minority atomic component defined by 
N^/Ntot- Fig. [T3] shows the contours of constant entropy assuming the system traverses 
through the adiabatic ramp, where the total entropy is given by S* = Sp{z^) + Sb{z^) + 

As discussed in the previous section, populations and adiabatic trajectories in 6-T plane 
assuming equal populations and different trap frequencies are plotted in Ref. p^] . Although 



we use different parameters from ones in Ref. 15|], characteristic behaviors do not change. 
Using this system, Ref. |l5| proposed cooling cycle with adiabatic sweep. 

The critical detuning 6c above which Bose condensation occurs at zero temperature is 
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FIG. 11: The phase diagram of the transition temperature of a bosonic minority component. We 
assume the ratio a = 3/4, and equal trap frequencies u)> = a;< = Wm- 



given by 



/ 1-a V^^ ( a V^^ 

^ 1 + a J \ 1 + a J 



At 5 — 5c and T = 0, all atoms in a minority component are converted into heteronuclear 
molecules since there is no bosonic thermal component 7V< = 0. At this point, populations of 
the majority component and the heteronuclear molecule are given by A^> = A^tot(l ~ ct)/(l + 
a) and = N^^^a/ (1 + The critical value 5c plays a role of matching the Fermi surfaces 
of a majority component and the heteronuclear molecule, i.e., = E^ + 5c, where and 
Ep are Fermi energies of the fermionic majority component and the fermionic heteronuclear 
molecule. With decreasing population of the minority bosonic component, the transition 
temperature of Bose-Einstein condensation Tc decreases. In the limit a — > 0, goes to 
zero, and the critical detuning 5c approaches the Fermi energy of the majority component 
5c Ep . Of course, the critical detuning 5c in the case a = has no physical meaning, 
because one has no minority bosonic atoms and no heteronuclear molecules. 

We note that the condition for a minority component BEG is given by /ii — —^2- Under 
this condition, the transition temperature Tc at zero detuning 5 = is given by 

ksTc = fuj^ -T7^-W7S7?^^tot . (48) 



,7^(l + 7;^)C(3)(l + a) 
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FIG. 12: (Color online) (A) The fraction of a fermionic majority component: Ny/Ntot- (B) The 
fraction of a bosonic minority component: A^</A'^tot- (C) The fraction of a fermionic heteronuclear 
molecular component: N^/Ntot- (D) The condensate fraction of a bosonic minority component: 
N^/Ntot. We assume the ratio a = 3/4, and equal trap frequencies (D> = u)< = Wm- 

From Eq. fHHj) . we find that the minority component does not undergo Bose condensation 
at zero detuning if the initial number ratio a is less than a critical value etc = ^>/{^> + ^'L)- 
We note that 6c is positive for a < ac- 

We first consider the case, a < where no condensation exists at 5 = 0. This means 
that all of the atomic minority component are converted into molecules at T = 0. In such 
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FIG. 13: (Color online) The contours of constant entropy. We assume the ratio a = 3/4, and 
equal trap frequencies a)> = a)< = Qjn- 



a case, the populations are given by 

iV>(5 = 0) = iV>,i„i-iV<, 



1 — a 
1 + a 



tot) 



iV<(5 = 0) = 0, 
N^{5 = 0) = 7V< 



a 



1 + a 



tot- 



(49) 



In this case, the maximum conversion efficiency Xo reaches 100% at T = 0. 

In the low temperature limit, the final entropy at zero detuning 5 = is given by 



= Sf{z>) + S'b(^<) + S'F(^m) 

-(iV>)2nV3 



6 



^B^^+4fcJ^) C(4) + ^B.^ 



1/3 



1 — a 



(iV„ 



1/3 



a 



\ 2/3 

^— iVtot + 7m ( ^tot 
1 + a ) \1 + a 



2/3' 



(50) 



where the approximation leading to the last line is valid when the final temperature satisfies 
the condition 



1 (1 
-n 



4C(4) \6 



1/6 
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-iVtot J(l-a)2/3 + Q;2/3 



1 + a J 



3/2 
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> ■ 



(51) 
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In a pure atomic gas before a sweep of the magnetic field, the initial populations are 
given by 



1 



>,ini , -""tot 

1 + a 



3l^>,mi;, 



a 



<,mi 



1 + a 



iVtot = iv; 



^B^ini 



(52) 



Qz{z<,\ 



The initial transition temperature Tc^ini and a condensate population iV^r^j^^ below Tc^mi are 
given by 



B-' c,mi 



a 



iV<,ini 

C(3) 



1/3 



1 + a 



tot 



^B^ini 



(53) 



C(3). 



In the low temperature limit, the initial entropy is given by 



S-f{z 



Fl,^>,iniJ 



S-niz^ 



^BTT^ 



{N: 



>,mi) 



'Bl,^<,ini 
\2 



/CB^ini , < , f ^B^ini 
+ 4fcB 



C(4) 



6 Vl + a 



tot 



1 1/3 



feB^ini 



(54) 



where the approximation leading to the last line is valid when the initial temperature satisfies 
the condition 



kBTini < 




(55) 



Connecting the final entropy at zero detuning 5 = with the initial entropy, we obtain 
the relation between the initial temperature Tjni and the final temperature Tf. In the low 
temperature limit, using Eqs. fl5Up and fl54p . we obtain the explicit formula 

1 



-T 



(56) 



;i -a)2/3 + a2/3- 
In the high temperature limit, the heteronuclear molecular population is so small that the 
contribution of the heteronuclear molecule entropy to the total entropy is small. As a result, 
one has Tf ^ Tjni in the high temperature region. Fig. [T3] (A) shows the numerical result for 
the final temperature as a function of the initial temperature, assuming a = 2/15 (which is 
less than ac = + lD^) = 1/2 for = ci}< = cjm). We also plot the result of the low 

temperature limit given in Eq. (!56|) . 
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We next consider the case a > ac = uj'^/iuj'^ + cuj^), where a condensate exists at zero 
detuning 5 = 0. The condensate population of the atomic minority component at zero 
detuning 5 = is given by 



;i + 7m)« 



^3 
/m 



C(3)7 



The population of each component below Tc at zero dutuning 5 = is given by 

1 1 



(57) 
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(58) 



iV„,(5 = 0) = iV<,i„i - iV< 



7^ 

/ m 



tot- 



1 + a 1 + 7i^ 

The population ratio of a fermionic majority component and a fermionic hetorunuclear 
molecule are constant with the bosonic minority component being Bose condensed, consistent 
with the fugacity relation z-^ = Zm at zero detuning 5 = below Tc where z^ = 1. 
The final entropy at zero detuning 5 = in the low temperature limit is given by 
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2A;R7r' ( - 



1/3 
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(59) 



where the approximation leading to the last line is valid for 
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:i+7 
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(60) 



By connecting the final entropy with the initial entropy, we obtain the relation between 
the initial temperature and the final temperature. In Fig.[Tl](B), we plot the numerical result 
for the final temperature assuming a = 3/4 (which is greater than ac = uj'^/{u!l,+uj'^) = 1/2 



for C(j> 
by 



ujra)- We also plot the analytical result in the low temperature limit, given 
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(61) 



Below Tc and at zero detuning 6 = 0, the population ratio of the majority component 
and the heteronuclear molecule becomes constant, which follows from the fugacity relation 
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FIG. 14: (Color online) The relation between the initial temperature and the final temperature 
at zero detuning 6 = 0. We assume equal trap frequencies In this case, a critical 

value Oc = 1/2 becomes a boundary, above which we have the condensation of the atomic minority 
component at zero detuning 5 = 0. Fig. (A) shows the relation with a = 2/15(< Oc = 1/2) and 
Fig. (B) shows the relation with a = 3/4(> Oc = 1/2). Solid lines are numerical results and dotted 
lines are results of low temperature limit. 

z> = Zai. As a result, the molecular conversion efficiency below Tc is given by 

^0 = - -3 ^ -3 - 62 

Fig. [15] shows the conversion efficiency as a function of the initial temperature. The 
solid line shows the conversion efficiency for a = 2/15, where Bose condensation does not 
occur. The dashed and dotted lines show the conversion efficiencies for a = 3/4 and 13/15, 
respectively, where a minority atomic gas becomes Bose-condensed at low temperature. A 
critical value ttc = ^^l/ ('^> + ^^m) (=1/2 for = lj^ = lD^) is a boundary that makes the 
trend of the conversion efficiency change, as shown in Fig. [151 

In Fig. [16], we plot the a-dependence of the maximum conversion efficiency for equal 
trap frequencies = lj^ = uj^. The maximum conversion efficiency Xo.max dramatically 
changes at a = 1/2. The maximum conversion efficiency Xo,max ranges from 50% to 100%, 
depending on the initial number ratio a. In general, for a > ac = LJy/{u't, + uj^), the 
maximum conversion is given by u"?^ / [a{u'^ + uj'^)] with a plateau. The maximum conversion 
efficiency ranges from Uy/{uj^ + o)^) to unity, as a function of the initial atomic fraction a. 
This result would not be expected if one applied the SPSS model as discussed in Sec. [TV] 
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Two cases {B> + F< Fm} and {F> + B< ^ Fm} are essentially the same. In the 
case {B> + F< ^ Fm}, maximum conversion efficiency with a plateau is given by Xo.max = 
ujp/{ujp + uj'^), where a trap frequency of the fermionic minority component uj^ is denoted 
as Up. On the other hand, in the case {F> + B< ^ F^}, maximum conversion efficiency 
for a > ac = COp/^iDp + cjj^) is given by Xo.max = ^p/l'^i^F + '^m)]) where a trap frequency 
of the fermionic majority component is denoted as up. We find that the maximum 
conversion efficiency is determined by initial number ratio and trap frequencies of fermionic 
components, which are both atoms and heteronuclear molecules. This result shows that 
fermionic components play a crucial role in determining each population, because energies 
of fermionic components are dominant for the total energy rather than bosonic one in the low 
temperature region. We note that one can achieve large conversion efficiency by making the 
trap frequencies of atoms larger than that of heteronuclear molecules, because the density of 
states of the heteronuclear molecule is higher than that of atoms in this case. This becomes 
more apparent when one writes the maximum conversion efficiency in Eq. (l62l) as 

7T- + o. A. , (63) 



where we define an atomic trap frequencies as cDa- This effect can be also seen in the case 
{B> + B< ^ Bin}, as well as cases {B> + F< F^} and {F> + B< ^ F^}. 



VI. GAS OF MAJORITY FERMIONIC ATOM, MINORITY FERMIONIC ATOM 
AND ITS HETERONUCLEAR MOLECULE 

Finally, we consider the case with a majority fermionic atom: {F> + F< ^ Bm}. Under 
the two constraints in Eqs. ([2]) and (jl]), the populations are given by 

^> = f^y^3(^>), iV< = (^X:F,iz^), = + (^^Qsiz^), (64) 

where is the population of the condensed heteronuclear molecule. The condition for 
heteronuclear molecular condensation is 2/ii + (1 — a)/i2 — 6 = 0. Fig. [T71 shows the phase 
diagram with the transition temperature where the heteronuclear molecules become Bose 
condensed. We assume a = 2/15, and equal trap frequencies lj^ = = u^a- Fig. [TS] shows 
fraction of each component. Figs. [18] (A), [18] (B) and [18] (C) show fractions of the majority 
atomic, the minority atomic and the molecular components defined by A^>/A^tot, -^</^tot 
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FIG. 15: (Color online) The molecular conversion efficiency as a function of tfie initial temperature. 
Wc assume equal trap frequencies In this critical value ac = 1/2 becomes 

a boundary, above which we have the condensation of the atomic minority component at zero 
detuning (5 = 0. We assume the ratio a = 2/15(solid line), a = 3/4(dashed line) and a = 
13/15(dotted line). 



o 

























































































■ >. - l_i I ■ - * - 1 1 





OM 0,S 0,75 



FIG. 16: The maximum molecular conversion efficiency Xo,max as a function of the initial popula- 
tion ratio of atoms a. We assume equal trap frequencies a;> = a;< = u)m- 
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FIG. 17: The phase diagram with the transition temperature of a bosonic minority component. 
We assume the ratio a = 2/15, and equal trap frequencies 

and N^/Ntot- Fig. [18] (D) shows the heteronuclear molecule condensate fraction defined by 
N^/Niot- Fig. [T9] shows the contours of constant entropy assuming the system traverses 
through the adiabatic ramp, where the total entropy is given hj S = S'f(-2>) + 5'f(-2<) + 

In Refs. ij] and [l6|, populations and adiabatic trajectories assuming population bal- 



anced gas and equal trap frequencies are plotted in 6-T plane. The characteristic behaviors 
plotted in this section do not change compared with results in Refs. [14] and [l^. 

The critical detuning Sc below which Bose condensation of molecules occurs at zero tem- 
perature T = is given by 

-I \ 1/3 / \ 1/3 



6, = nLo^[6-—Ntot] +to< 6^— iVtot . (65) 

At 6 = Sc, heteronuclear molecules do not exist since there is no thermal component at 
T = 0. At this point, populations of the majority and the minority components are given 
by A^> = A^tot/(l + «) and A^< = A^tot«/(l + «)• The critical value is also represented as 
6c = Ep +Ep, where Ep and Ep are Fermi energies of a majority and minority components. 
Just below 5c at T = 0, although the single-particle lowest energy of heteronuclear molecules 
is higher than that of atoms, condensed heteronuclear molecules exist owing to Pauli's 
exclusion principle. 
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FIG. 18: (Color online) (A) The fraction of a fermionic majority component: Ny/Ntot- (B) The 
fraction of a fermionic minority component: N^/Ntot- (C) The fraction of a bosonic heteronuclear 
molecular component: N^/Ntot- (D) The condensate fraction of a bosonic heteronuclear molecular 
component: N^/Ntot. We assume the ratio a = 2/15, and equal trap frequencies 

At zero detuning 5 = 0, the condition for heteronuclear molecular condensation is given 
by /ii = —(1 — a)/i2/2. Chemical potentials of majority and minority components /i> and 
/i< are given by /i> = — (1 + a) ^2 and /i< = (1 + a)fi2- At T = 0, a part of atoms in a 
majority component is not converted into heteronuclear molecules, which forms the Fermi 
surface. From this condition, the chemical potential of a majority component should be 
positive, i.e. /i2 < 0, and thus one has no minority component due to the negative chemical 
potential at T = and 6 = 0. The populations at T = and zero detuning 6 = are given 
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FIG. 19: (Color online) The contours of constant entropy. We assume the ratio a = 2/15, and 
equal trap frequencies a)> = a;< = Wm- 
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At T — 0, the conversion efficiency xo is 100%, which gives the maximum conversion. 
The final entropy at zero detuning 5 = in the low temperature limit is given by 
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where the approximation leading to the last line is valid for 
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On the other hand, the initial entropy in the low temperature limit is given by 

'S'ini = >S'> ini + *S'< ini 
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(69) 



Connecting the final entropy at zero detuning with the initial entropy, we obtain the 
relation between the initial temperature Tini and the final temperature Tf. In the low tem- 
perature limit, the analytical expression for this relation is given from Eqs. (!67|) and (!69|) 

as 



l_a)2/3^-i- 



(70) 



In the high temperature limit, the heteronuclear molecular population is so small that the 
contribution of the heteronuclear molecule entropy to the total entropy is small. As a result, 
one has Tf ^ Tini in the high temperature region. Fig. [20] shows the numerical result for 
the final temperature as a function of the initial temperature for a = 2/15. We also show 
the analytical result in the low temperature limit given in Eq. (170]) . In the population 
balanced case a = 1, the final entropy in the low temperature limit is only composed of 
the heteronuclear molecule entropy given by Sf = 4A;B7m (^BTf/^>)^ C(4)- In this case, the 
relation between the initial and final temperatures in the low temperature limit is given by 



fceTf = - — 
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1/3 
tot 
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(71) 



2 7n, LC(4)3i 

Fig. [21] shows the conversion efficiency as a function of the initial temperature given 
from the heteronuclear molecular population at zero detuning 6 = 0. One can see that the 
maximum conversion reaches 100% at T = 0. Molecular conversion efficiencies in trapped 
gases composed of two component Fermi atoms assuming equal trap frequencies and equal 
populations are discussed in Refs. Hand [l4l |. The behavior shown in Fig. [21] is qualitatively 



consistent with the results in Ref. 



VII. EXPLICIT FORMULA FOR CONVERSION EFFICIENCY 

In Ref. [it, the SPSS model predicted the conversion efficien^ as a function of the 
peak phase space density. On the other hand, Williams et al. [Tj explicitly derived the 
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FIG. 20: (Color online) The relation between the initial temperature and the final temperature 
at zero detuning 6 = 0. The solid line shows the result of the numerical calculation. The dotted 
line shows the result of the low temperature limit. We assume the ratio a = 2/15, and equal trap 
frequencies a)> = a;< = Wm- 




FIG. 21: The molecular conversion efficiency as a function of the initial temperature. We assume 
the ratio a = 2/15, and equal trap frequencies a;> = a)< = Qui- 
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conversion efficiency as a function of the initial peak phase space density within the classical 
gas approximation. Both results are in good agreement with the experimental result In 
this section, we explicitly derive the conversion efficiencies for the heteronuclear Feshbach 
molecule in population imbalanced CclSG clS db function of initial peak phase space density 
above T^, treating a majority component as a quantum degenerate gas. 

We use the Bose or Fermi distribution function for a majority component, while we use 
the Maxwell-Boltzmann distributions for a minority atomic component and heteronuclear 
molecules, because these minority densities are so small that interparticle distances are larger 
than the thermal de Broglie wavelength. In this approximation, the Bose and Fermi integrals 
for the minority component and heteronuclear molecules are given approximately by 



z. 



(72) 



First, we consider the case with majority bosonic atoms within the approximation men- 
tioned above, the number of each component is given by 
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(73) 



The entropy of each component is given by 
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S< = A;BiV<(4-lnz<), 
Sm = kBN^ (4: - In z^) . 

The molecular conversion efficiency is then given by 
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where t = k-^T / {fvjjyN^ 
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According to the above equation and two constraints in Eqs. ([2]) and (jlj) for populations, 
we obtain Zm, 2;<, and t^Qz^Zy) at zero detuning 5 = as a function of xo; 
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Prom the relation Zm — z>z<^ for fugacities at zero detuning 5 = 0, the fugacity of the 
majority component 2;> as a function of xo is given by 



Xo 



z< 1 - Xo V7,'^ 



(77) 



Using these equations, we obtain the final total entropy at zero detuning 5 = as a function 
of a, Xo, 7< and 7^, 
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On the other hand, the initial numbers of particles are given by 
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(79) 



A^<,ini = 
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tot- 



The peak phase space density of a bosonic majority component is given by 

Ppk = Atn>(r = 0) = e;3/2(^>), 



(80) 



where A> is a thermal de Broglie wavelength of a majority component given by A> = 
[27rh^ / (m^kBT)]^^'^ , with m> being an atomic mass of the majority component. We obtain 
the fugacity z> as a function of p^^; 



(81) 



By making use of the equation of the second constraint = -/V<,ini — Q;A^>,ini, we obtain 
the relation between 2;<,ini and 2;>,ini; 



7<^<,ini = Q;^3(^>,ini)- 



(82) 
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The initial total entropy is given as a function of the ratio a and the initial peak phase 
space density of a majority component Ppkini! 

Q. . — Q . . _|_ Q . . 
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(83) 



By equating the final entropy with the initial entropy, we obtain the relation between the 
molecular conversion efficiency xo and the initial peak phase space density of a majority 
component Pp^in; as 
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Since the fugacity of the majority component z> should be less than unity, it follows from 
Eq. (!77|) that the range of the conversion efficiency above the transition temperature of the 
majority atomic component BEC is limited as 



,-,3 



XO < XO,max = -o ■ (85) 

This maximum conversion efficiency is consistent with our results discussed in Sec lIIII andllV t 
in the case where the majority component is bosonic. 

Using the similar procedure, we consider the case with a majority fermionic component. 
In this peak phase space density is given by 

Ppk,ini = Air2>(r = 0) = ^3/2(^>). (86) 

The relation analogous to Eq. (18^ for a fermionic majority component is obtained by simply 
replacing the Bose integral with the Fermi integral. 
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According to Eqs. (18^ and (IHTl) . we find that the initial atomic population ratio a 
and the ratio of the trap frequencies of the minority atomic component and hetoronuclear 
molecular component uo-^Juo^ as well as the initial peak phase space density of the majority 
atomic component Pp^jn; determine the molecular conversion efficiency Xq iii the adiabatic 
ramp to form the heteronuclear Feshbach molecule. 

In the high temperature region, one can also use the Maxwell-Boltzmann distribution 
for a majority component; Qn{,z>) ~ -2> or J-'n{z^) ~ In this situation, a fugacity of a 
majority component is approximately given by 

^> ~ Ppk,im- (88) 

Quantum statistics is not important in this region. For classical gases, the relation between 
the molecular conversion efficiency ^^nd the initial peak phase space density of a majority 
component p^^^ j^^- is given by 



^ Xo + In — - = lnp> (89) 

If we consider the population balanced case a = 1 and equal trap frequencies cc;< = ujm, we 
recover the result given by Williams a/ [3]; 

lnp>,i„i = 2xo + ln(^^). (90) 

By using the relation 7<Ppkini = "^Ppkini obtained from Eq. (jlj) for the initial state, we also 
obtain the relation between the molecular conversion efficiency xo and the initial peak phase 
space density of a minority component Ppy-^^i', 



1 / - \ 3 



mi 



(91) 
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Fig. [22] shows the molecular conversion efficiency Xo a function of the initial peak 
phase space density Pp^ini- We assume a = 2/15, and equal trap frequencies 
Fig. [22] (A) shows the case where the majority atomic component is bosonic. The solid line 
is the numerical result of{B> + B<^ B^} discussed in Sec. IIIII The dashed line is the 
numerical result of{B> + F<^ Fm} discussed in Sec. [IV] The dot-dashed line is the result 
of our formula given in Eq. ([84|) . reproducing two numerical lines quite well. The dotted 
line is the result of the classical limit given in Eq. ([891) . which overlaps with three lines in 
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FIG. 22: (Color online) The molecular conversion efficiency as a function of the initial peak 
phase space density of a majority component. We use the parameter a = 2/15, and equal trap 
frequencies. (A) The solid and the dashed lines are the numerical results of {B> + B< ^ Bm} 
and {B> + F < <— > Fm}- The dot-dashed line is the result of our formula given in Eq. (j84p . Those 
lines agree quite well with each other. The dotted line is the result of the classical limit given in 
Eq. ([89]) . (B) The solid and the dashed lines are the numerical results of {F> + B< <-5- Fm} and 
{F> + F< ^ Bm}- The dot-dashed line is the result of our formula given in Eq. (|87p . Those lines 
agree quite well with each other. The dotted line is the result of the classical limit given in Eq. 

(El. 

the low initial peak phase space density limit. Fig. [221 (B) shows the case where the majority 
atomic component is fermionic. The solid line is the numerical result of {F> + B< ^ Fm} 
discussed in Sec. [V| The dashed line is the numerical result of {F> -|- F< ^ Bm} discussed 
in Sec. I VII The dot-dashed line is the result of our formula given in Eq. ([87]) . reproducing 
two numerical lines quite well. The dotted line is the result of the classical limit given in 
Eq. ( [89l) . which agrees with three lines in the low initial peak phase space density limit. 

Only in the case {F> + B< <-> Fm} discussed in Sec. [Vj the maximum conversion efficiency 
depends on the initial atomic population ratio a. Fig. [251 shows the conversion efficiency as 
a function of the initial peak phase space density of the fermionic majority component in the 
case {F> + B< Fm}, assuming equal trap frequencies = c<}< = cUm and a = 3/4 > ac(= 
1/2). In this case, the maximum conversion efficiency is given by l/(2a). The solid line is 
the numerical result discussed in Sec. [V| The dashed line is the result of our formula given 
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in Eq. (187|) . The dotted line is the result of the classical limit given in Eq. (l89l) . According 
to Fig. [231 our formula of Eq. (187]) does not reproduce the numerical result in the high initial 
peak phase space density region, differing from the case a = 2/15(< = 1/2) in Fig. [22] 
(B). In the step deriving Eq. fISTI) . we treated that the minority atomic component and 
molecules as classical Maxwell-Boltzmann gases. For a > ac (= 1/2), this approximation 
breaks down in the region where the quantum degeneracy of the minority atomic component 
cannot be ignored. The emergence of Rose statistics that the minority atomic component 
resides causes the disagreements between the numerical results and our formula as shown in 
Fig. [231 

To summarize, we found that the molecular conversion efficiency in an adiabatic sweep is 
well described as a function of the initial number ratio and the ratio of the trap frequencies 
of the minority atomic component and the heteronuclear molecular component as well as 
the initial peak phase space density 

XO = XO ( Ppk,ini. ^ ) ' (92) 

in the case where populations of the minority component and the heteronuclear molecule 
are small. 



VIII. CONCLUSION 



In this paper, we studied formations of heteronuclear Feshbach molecules in population 

n 

imbalanced atomic gases, extending the recent work [7| on the Feshbach molecule formation. 
At low temperature in quantum degenerate regime, quantum statistics of atoms plays an 
important role in determining conversion efficiencies. 

When the majority and minority atomic components are both bosonic, the maximum con- 
version efficiency is determined by the trap frequencies of the minority atomic component 
and heteronuclear molecules. Our calculation is in good agreement with the recent experi- 
ment Sj without any fitting parameters. An important finding is that one cannot convert 
any atoms into heteronuclear Feshbach molecules in the limit Tin; — 0. On the other hand, 
when gases are composed of fermionic atoms and bosonic atoms, the maximum conversion 
efficiencies are determined by the trap frequencies of fermionic atoms and heteronuclear 
molecules as well as initial number ratio. In the case that both atomic components are 
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FIG. 23: The molecular conversion efficiency as a function of the initial peak phase space density 
of a majority component, in the case {F> + B< <-^> Fm} assuming a > Oc- We assume a = 3/4(> 
Oc = 1/2) and equal trap frequencies. The solid line is the numerical result. The dashed line is the 
result of our formula given in Eq. (j87p . The dotted line is the result of the classical limit given in 
Eq. dMI)- 

fermionic, the maximum conversion efficiency is 100%. In general, when atoms are not Bose 
condensed even at zero temperature at zero detuning, the molecular conversion efficiency 
reaches 100%. When one atomic component undergoes Bose-Einstain condensation but the 
other component does not at zero detuning, the conversion efficiency does not reach 100%, 
exhibiting a plateau. 

In the region where the gases are not condensed, the conversion efficiency is described 
as an explicit function of the initial number ratio of atoms and trap frequencies of minority 
component and heteronuclear molecule as well as the initial peak phase space density of a 
majority atomic component. We found that in the low-density region where Bose-Einstein 
condensation does not appear, the conversion efficiency is a monotonic function of the peak 
phase space density, but independent of statistics of the minority component. 

Throughout this paper, we assumed equal trap frequencies in all figures for simplicity. 
This simple assumption is valid for loosely bounded Feshbach molecules composed of isotopic 
atoms. Although qualitative behaviors with respect to all figures do not change when we 
assume different trap frequencies, we note that the maximum conversion efficiencies depend 
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on trap frequencies if BEC of the single atomic component appears. 

'. ■i'inally, we note that the theory of Feshbach molecule formationproposed by Williams et 
al. 3] brings the result quite different from the SPSS model of Ref. jl] in the low temperature 
region. 
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